Generating Many Majorana Modes via Periodic Driving: A Superconductor Model 
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Realizing Majorana modes in condensed-matter systems is of vast experimental and theoretical 
interests, and some signatures of Majorana modes have been measured already. To facilitate future 
experimental observations and to explore further applications of Majorana modes, generating many 
Majorana modes at ease in an experimentally accessible manner has become one important issue. 
This task is achieved here in a one-dimensional p-wave superconductor system with the nearest- and 
next-nearest-neighbor interactions. In particular, a periodic modulation of some system parame- 
ters can induce an effective long-range interaction (as suggested by the Baker-Campbell-Hausdorff 
formula) and may recover time-reversal symmetry already broken in undriven cases. By exploiting 
these two independent mechanisms at once we have established a general method in generating 
many Floquet Majorana modes via periodic driving. 

PACS numbers: 03.67.Lx, 03.65. Vf, 71.10.Pm 



Introduction. — The Majorana fcrmion, a particle of its 
own anti-particle [1], is attracting tremendous attention 
[2-5]. In addition to its fundamental interest [6-9], its 
potential applications in topological quantum computa- 
tion are also noteworthy [10, 11]. Along with consider- 
able theoretical studies [7, 10, 12-20], the experimental 
search for Majorana modes in condensed-matter systems 
has become a timely and important research topic. In- 
deed, following the theoretical results in Refs. [15, 16, 21- 
25], the zero-bias conductance peaks observed recently 
[26-29] are regarded as a signature of Majorana modes 
in one-dimensional (ID) spin-orbit coupled semiconduc- 
tor nanowires. However, the observed zero-bias peaks can 
be due to other reasons as well, e.g., the strong disorder 
in the nanowire [30, 31] or smooth confinement potential 
at the wire end [32]. This being the case, the formation 
of Majorana modes in these nanowire systems are yet to 
be double-confirmed by other approaches. 

To facilitate experimental observation of Majorana 
modes, their signal strength should be enhanced so as 
to combat against experimental disorder [24, 30, 33] and 
contaminations from thermal excitations [21-25]. To that 
end, it is constructive to find a general method to form 
many Majorana modes within one single system [34-36]. 
This turns out to be a challenging task. One main rea- 
son is that the formation of many Majorana modes needs 
the protection of time-reversal symmetry [36-38], thus 
imposing a strong restriction. Reference [37] remarked 
that three or even more pairs of Majorana modes can be 
obtained by considering a "longer-range" interaction in 
a system. This is a stimulating suggestion, hinting that 
the more long-ranged the interaction is, the more Ma- 
jorana modes one might obtain. Generating many Ma- 
jorana modes then reduces to the following theoretical 



question: how to synthesize a long-range interaction in 
a topologically nontrivial condensed-matter system while 
maintaining time-reversal symmetry? 

Our answer to the above question is rather straight- 
forward. Given that periodic driving is becoming one 
highly controllable and versatile tool in generating dif- 
ferent topological states of matter [39-44], we advocate 
to use a periodic driving protocol to induce effective long- 
range interactions and to restore time-reversal symmetry 
(if it is broken without driving). Specifically, we switch 
a Hamiltonian from Hi for the first half period to H 2 for 
the second half period. The Floquet operator U, i.e., the 
unitary evolution operator associated with each modula- 
tion period T, is then given by 



U(T) = 



(1) 



where an effective Hamiltonian H c fi for the modulated 
system is also defined. Using the Baker-Campbell- 
Hausdorff (BCH) formula, one finds that H c s defined in 
Eq. (1) is formally given by 



Hi H 2 
T 2 



iT 
8h 



[H 2 ,Hi 



96/T 



[{H 2 - Hi),[H 2 ,Hi] 



(2) 



Clearly then, even if Hi or H 2 are short-range Hamilto- 
nians, the engineered H c g may still have long-range hop- 
ping or pairing terms via the nested-commutator terms 
[see Eq. (2)]. Thus, the remaining job is to design such a 
protocol so that H e g also possesses time-reversal symme- 
try. Interestingly, in the first proposal to realize Floquet 
Majorana modes [44], no more than two pairs of Ma- 
jorana modes can be generated precisely because time- 



2 



reversal symmetry is not restored by the periodic driving 
therein. 

In the following we present our detailed results using 
a model of a ID spinless p-wave superconductor with the 
nearest- and next-nearest-neighbor (NNN) interactions 
only. The introduction of the NNN interaction will fa- 
cilitate us to have time-reversal-symmetry preserved and 
broken cases and thus to realize their conversion by our 
proposed driven protocol. Under open boundary condi- 
tion and under a periodic modulation of superconducting 
phase parameters, we not only demonstrate that many 
Floquet Majorana modes (e.g., 13 pairs in one case) can 
be generated, but also show that the exact number of the 
Majorana modes may be widely tuned by scanning the 
modulation period. These results also shed more light on 
the inherent advantages of driven systems in synthesizing 
Majorana modes. 

Static model. — We start from the Kitaev model Hamil- 
tonian for a ID spinless p-wave superconductor 

N N-l 

H = -nY^ c \ c i- X](* lC i c '+ 1 + AlC i c m +h - c -) 
;=i i=i 

- ^4 c i+2 + a 2 44 +2 + h.c), (3) 
i=i 

where [i is the chemical potential, t a and A a = |A a | e % ^ a 
with a = 1 (a = 2) describe the nearest- (next-nearest-) 
neighbor hopping amplitude and pairing potential respec- 
tively, and 4> a is the associated superconducting phases. 
All energy- related parameters are scaled by |Ai| and 
hence dimensionless. For convenience we set ft = 1 in 
our detailed calculations. Majorana operators here re- 
fer to (c; + cj) or i(ci — cj). Such synthesized Majorana 
modes may appear as edge modes under open bound- 
ary condition, if the bulk band structure is topologically 
nontrivial. 

The relative phase <j) = 0i — 02 determines the topolog- 
ical class of H [45]. For = and n, H has time-reversal 
and particle-hole symmetries. These cases then belong to 
the so-called "BDI" class characterized by a topological 
invariant Z. For other values of 0, H has particle-hole 
symmetry only and falls into the so-called "D" class char- 
acterized by a topological invariant Zi- The D class can 
generate at most one pair of Majorana modes. As to the 
BDI class, despite its potential in forming many Majo- 
rana modes [38], at most two pairs of Majorana modes 
can be generated here due to the short-range nature of 
if. 

Driven model. — We now turn to periodically driven 
cases under a driving protocol given by Eq. (1). The 
emergence of Floquet Majorana modes is directly con- 
nected with topological properties of the eigenstates of 
the Floquet operator U{T). Let \u) be an eigenstate 
of U(T) with an eigenvalue e~ teT , namely U(T)\u) — 
e~ leT \u). Evidently, this eigenvalue index e is defined 




FIG. 1: (Color online) (a) Winding of the unit vector n(k) 
[see Eqs. (5) and (6)] as k varies from — -n to n, in different 
parameter regimes under our driving protocol. t\ = 1; t% = 5 
for W = -3, h = 3 for W = -2, t 2 = for W = 0, and 
t2 = —3 for W = 1. (b) Topological phase diagram of our 
driven system plotted in the (ti,fo) plane. Indicated on each 
regime is the winding number W, which gives the number of 
pairs of Majorana modes. The solid and dotted line indicates 
a gap closing (of Ek) at k — or ±7T, while the dashed line 
corresponds to a gap closing at k = n/2. Other parameters 
are fi = -10, |A 2 | = 2.5 and T = 0.2. 

only up to a period 2ir/T and hence called "quasi- 
energy" . The periodicity in e may lead to a novel topo- 
logical structure in driven systems, with the correspond- 
ing topological classification revealed by the homotopy 
groups [42]. However, if the driven system belongs to a 
trivial class associated with this novel topological struc- 
ture, then topological properties of the driven system is 
fully characterized by H e g defined in Eq. (1) [38]. This 
will be the case for our driving protocol proposed below. 

Specifically, we switch between the two Hamiltonians 
Hi and H2 by the following: in the first half period, Hi = 
H(4>i,4>2) with both superconducting phase parameters 
4>\ and (f>2 fixed; whereas in the second half period, we 
swap 0i and <p2 so that H2 = H(<f>2,<f>i). Without loss 
of generality we choose <\> = <fii — <f>2 = ±n/2. Thus, 
within each half period, the Hamiltonian is in class D 
that breaks time-reversal symmetry. In addition to a 
possible generation of long-range interactions for H e g, 
this driving protocol is designed to recover time-reversal 
symmetry. In particular, let K. be a conventional time- 
reversal operator and G = e~ l 2 ^i c i Cl be a gauge 
transformation operator. Considering a generalized time- 
reversal operator K = KG, we find 

ICUiT)^- 1 = e^e^ = U\T). (4) 

This constitutes a direct proof that our driven system 
now possesses time-reversal symmetry, and as a result 
its topological class is switched from class D to class 
BDI. To further examine this restored time-reversal sym- 
metry, we work in the momentum representation and 
directly find an analytical H c g from Eq. (1). We de- 
fine Cfe = ^}2 l cie~ lkl / y/N and introduce the particle- 
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hole representation C k — [cfc,c_ fe ] T . A standard pro- 
cedure then leads to H e g = X^^gbz ClH a ff(k)C k , with 
H c g(k) = E k ft(k) ■ a, where a represents the Pauli ma- 
trices [41]. The three components of n(k) are given by 
ni(fc) = 0, and 



(a) 



(b) 



n 2 (k) 
n 3 (k) 



9i,k sin(sfcT) _ 2g 2: kVk sin 2 s k T/2 
s 2 k sm(E k T) ' 
2ffi,fe52,feSin 2 s k T/2 



s k sin(E k T) 
rjk sin(sfcT) 
s k sin{E k T) 



; sin(£ , fe T) 



(5) 
(6) 



where g a , k = |A a |sin(afc), s k = (vt + Y, a 9l, k ) 1/2 , 
Vk = — M - 2J2a t aCOs(ak), and cos(E k T) = cos(s k T) + 
2{g^ k / s\)s\v? s k T /2. For each value of k, one obtains 
two values of E k and hence two values for the quasi- 
encrgy e. Consistent with the K, symmetry, we now have 
H* s (—k) = H c s(k). Noting the inherent particle- hole 
symmetry of i? c ff, one may construct a chiral symmetry 
for H c ff , a fact consistent with our above result that n(k) 
is in the yz plane for all k. The above analysis makes it 
clear that our driving protocol changes both the under- 
lying symmetry and the topological class of the system. 

Without a gap closing between the two branches of E k , 
the topological invariant Z in class BDI can be obtained 
by the integer winding number W = -s-k € Z, where 
k = arctan[n3(fc)/rj2(fc)]. Geometrically, W means the 
number of times the vector n(k) rotates around the ori- 
gin point as k varies from — 7r to n. An example is shown 
in Fig. 1(a). The number of pairs of Majorana modes 
under open boundary condition is then given by \W\ . As 
some system parameters continuously change, gap clos- 
ing and consequently topological phase transitions occur 
[7]. Figure 1(b) depicts a phase diagram, obtained by 
explicitly evaluating W. It is seen that that \W\ ranges 
from to 3. This indicates that three pairs of Majorana 
modes can be formed in our driven system. This is be- 
yond the expectation for the undriven model, where the 
NNN interaction can give at most two pairs of Majorana 
modes. Therefore, the finding of \ W\ = 3 in some param- 
eter regime is the first clear sign that our driving protocol 
may synthesize some features absent in the static model. 
The boundaries between different topological phases of 
our driven system, as presented in Fig. 1(b) by the solid, 
dashed, and dotted lines, are also interesting on their 
own right. For example, the solid and dotted lines de- 
pict the topological phase transition points at which W 
jumps by one. This is found to go with the gap clos- 
ing at k — or ±7r. Noting that n 2 {k) = for k = 
or ±7r, the solid and dotted lines are found to satisfy 
sin[(±2ti - 2t 2 - ^)T] = 0, which gives n 3 (k) = for 
k = or ±7r. On the other hand, the dashed line in 
Fig. 1(b) gives the phase transition points at which W 
jumps by two. This happens at k = n/2 and the critical 
t 2 values satisfy t 2 = (^+ v /(tt/T) 2 - |Ai| 2 )/2, for which 
n 2 (k) = n 3 {k) = for k = tt/2. Such a phase border line 
is independent of t% and IA2I, thus explaining why the 
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FIG. 2: (Color online) Comparison between quasi-energy 
spectrum for the driven case (a) and energy spectrum for the 
static case (b), both numerically obtained under open bound- 
ary condition. The dashed (blue) line in (a) stands for two 
degenerate pairs of Majorana modes, and the solid (black) 
lines in (a) and (b) stand for a single pair. The results agree 
with the topological phase diagram in Fig. 1(b). Here ti = 1, 
N = 200, and other parameters are the same as in Fig. 1(b). 



dashed line in Fig. 1(b) is horizonal. 

To confirm our theoretical results presented in Fig. 1 
we carry out numerical calculations of the quasi-energy 
spectrum e under open boundary condition. Because 
e = tt/T is equivalent to e = —ir/T, Floquet Majorana 
modes have two flavors: one at e = and the other at 
e = ±ir/T. The second flavor is certainly absent in an 
undriven system [44]. For fixed ti = 1 and a varying t 2 , 
Fig. 2(a) depicts the formation of both flavors of Floquet 
Majorana modes obtained numerically, with the second 
flavor emerging in a wider parameter regime. The total 
number of pairs of Majorana modes should equal \W\ (if 
the winding number is well defined). As an example, con- 
sider the case with t\ = 1 and t 2 = 4. Figure 2(a) shows 
the formation of two degenerate pairs of Majorana modes 
at e = ±ir/T and one pair of Majorana mode at e = 0. 
This agrees with the W = —3 region shown in Fig. 1(b). 
Likewise, all other details shown in Fig. 2(a) are fully con- 
sistent with our analytical results shown in Fig. 1(b). We 
have also studied the dynamics of the formed Majorana 
modes in one full period of driving: they are indeed well 
localized at two edges. Further, as a comparison with 
our static model H, wc plot in Fig. 2(b) our system's en- 
ergy spectrum in the absence of driving. It is seen that at 
most one pair of Majorana modes can be formed. Indeed, 
theoretical analysis for the undriven case in the class D 
yields that the associated Z 2 topological number is given 
by v = sgn[(2t2+M) 2 — 4ii] and at most one pair of Majo- 
rana modes may be formed when v = — 1 . For a negative 
/i but with large \p\ , one has v = — 1 only in a very narrow 
parameter regime. The parallel driven case with a large 
large is however different: one may still obtain three 
pairs of Majorana modes. Thus, even in regimes of large 
I /it I, our driving protocol can still generate more Majo- 
rana modes than the static case. This is both interesting 
and useful because in general, the regime of large is 
preferred for the protection of Majorana modes against 
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(a) (b) 




FIG. 3: (Color online) The expansion coefficients of H e g for 
T = 0.2 (a) and T = 2.0 (b), and of the static H for the real 
(c) and imaginary (d) parts, when H e fi or H are expanded as 
a quadratic function of the operators (ci, • • • ,cjv, &[,••• ,cjy), 
with N = 200, and I and m being the operator base indices 
ranging from 1 to 400. Other system parameters are the same 
as in Fig. 2. 



\fi\, in Fig. 3(c). This further explains why a driven sys- 
tem may generate many Majorana modes despite a large 
\fx\ in the undriven model. 

Results in Fig. 3 motivate us to explore the formation 
of Floquet Majorana modes with sufficiently large values 
of T. There is also one twist as we increase T. That is, 
the quasi-energy gap may be generically closed at e = 
but still open at e = ±tt/T. Consequently the winding 
number W is no longer well-defined and one is forced to 
count the formed pairs of Majorana modes, one by one, 
at e = ±7r/T. Wc indeed do this and present in Tab. I 
the number of pairs of Majorana modes we obtain, for an 
increasing T and for different choices of £2. For T = 0.5, 
the best observation is the generation of 4 pairs of Majo- 
rana modes. For T = 1.0, it is possible to achieve 7 pairs. 
For T = 2.0, as many as 13 pairs of Majorana modes can 
be formed. Interestingly, as indicated in Tab. I, for large 
values of T such as T = 1.0 and T = 2.0, our driving 
protocol can also form several pairs of Majorana modes 
with £2 = 0. Note also that as more Majorana modes 
are generated by increasing T, the bulk quasi-energy gap 
at s = ±7r/T decreases in general, leading to a larger 
"penetration length" (into the bulk) for the synthesized 
Majorana modes. Considering the necessary protection 
of Majorana modes by a nonzero bulk gap, one may not 
wish to push our driving protocol too far. 



strong disorder in actual experiments. 

In efforts to generate even more Majorana modes, we 
now extend our direct numerical studies to other param- 
eter regimes. Remarkably, the BCH formula in Eq. (2) 
indicates that as T increases, the nested commutators on 
the right hand side of Eq. (2) will have heavier weights. 
An increasing T can then induce longer-range interac- 
tions in H B g. This trend is investigated in Fig. 3, where 
the expansion coefficients of if e ff [numerically obtained 
from Eq. (1)], with H e g expanded as a quadratic func- 
tion of the operators (ci, • • • , Cn, c\, ■ ■ ■ , cjy) T , are shown 
for two different values of T. For comparison, the expan- 
sion coefficients for the static case are also plotted in 
Fig. 3(c,d). A few interesting observations can be made 
from Fig. 3. First, the plotted expansion coefficients of 
fleff are all real, which is different from the shown static 
case with both real and imaginary coefficients. This dif- 
ference reflects the restored time-reversal symmetry for 
the driven case. Second, in sharp contrast to the results 
shown in Fig. 3(c,d), coefficients for quite long-range 
hopping/pairing (e.g., across more than 10 sites) can be 
appreciably nonzero for iJ g in both cases of T = 0.2 
and T = 2.0. The latter case, plotted in Fig. 3(b) with 
wider stripes, confirms the emergence of longer-range 
terms with considerable weights as T increases. Third, 
the diagonal terms in the expansion shown in Fig. 3(a,b) 
(which can be understood as an effective chemical poten- 
tial) are much smaller than the diagonal elements, i.e. 



TABLE I: Number of Majorana modes localized at each 
boundary for different values of T. Other parameters are the 
same as in Fig. 2(a). The parenthese indicate that the system 
is in a close neighborhood of a topological phase transition. 
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Conclusions. — A periodic driving has the capacity to 
restore time-reversal symmetry and to induce an ef- 
fective long-range interaction. With these two mecha- 
nisms working at once, the generation of many Majorana 
modes is achieved using a standard p-wave superconduc- 
tor model under certain periodic modulation. 

In terms of possible experimental confirmation of our 
predictions, our model may be realized with cold atoms 
or molecules in a designed optical lattice. Explicitly, the 
nearest and NNN hopping (£„) can be realized by a simple 
zigzag chain lattice [35], with the hopping strength ad- 
justable by the lattice geometry. The chemical potential 
(/i) is controllable through the optical trap potential or a 
radio frequency detuning. The pairing terms (A Q ) may 
be induced by a Raman induced dissociation of Cooper 
pairs forming an atomic BCS reservoir and the associ- 
ated superconducting phases (<f> a ) can be tuned by com- 
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plex Rabi frequencies [44]. An alternative experimental 
realization is to use the recently proposed quantum-dot- 
superconductor arrays with a zigzag geometry [46] . Here 
the chemical potential can be gate controlled. The pair- 
ing potentials can be proximity-induced and their phases 
can be tuned via applying fluxes on the superconducting 
islands. The Majorana modes formed in our system may 
be probed using techniques analogous to what is being 
used for undriven systems [35], but now with the hope 
of some enhanced signals if a measurement exploits the 
simultaneous generation of many Majorana modes. 
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